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ABSTRACT 

A characteristic feature of the Fourier Power Density Spectrum (PDS) ob- 
served from black hole X-ray binaries in low/hard and intermediate spectral 
states is a broad band-limited noise, characterized by a constant below some fre- 
quency (a “break” frequency) and a power law above this frequency. It has been 
shown that the variability of this type can be produced by the inward diffusion 
of the local driving perturbations in a bounded configuration (accretion disk or 
corona). In the framework of this model, the perturbation diffusion time to is 
related to the phenomenological break frequency, while the PDS power-law slope 
above the “break” is determined by the viscosity distribution over the configura- 
tion. The perturbation diffusion scenario explains the decay of the power of X-ray 
variability observed in a number of compact sources (containing black hole and 
neutron star) during an evolution of theses sources from low/hard to high/soft 
states. We compare the model predictions with the subset of data from Cyg X-l 
collected by the Rossi X-ray Time Explorer (RXTE). Our extensive analysis of 
the Cyg X-l PDSs demonstrates that the observed integrated power P x decreases 
approximately as a square root of the characteristic frequency of the driving os- 
cillations v dr . The RXTE observations of Cyg X-l allow us to infer P dr and t 0 as 
a function of u dr . Using the inferred dependences of the integrated power of the 
driving oscillations P dr and to on u dr we demonstrate that the power predicted 
by the model also decays as P x ,diff oc z /^. 0 ' 5 that is similar to the observed P x 
behavior. We also apply the basic parameters of observed PDSs, power-law index 
and low frequency quasiperiodic oscillations, to infer Reynolds (Re) number from 
the observations using the method developed in our previous paper. Our analysis 
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shows that Re— number increases from values about 10 in low/hard state to that 
about 70 during the high/soft state. 

Subject headings: accretion, accretion disks — black hole physics — stars:individual 
(Cyg X-l) :radiation mechanisms: nonthermal — physical data and processes 


1. Introduction 

The main goal of the presented work is to explain a decay of the emergent time variability 
of X-ray emission in compact sources when these sources evolve from low/hard (LH) to 
high/soft (HS) states [see Remillard & McClintock (2006) and Titarchuk, Shaposhnikov & 
Arefiev (2007), hereafter TSA07, for details of observations]. In particular, power density 
spectrum (PDS) of black hole binaries in hard states is dominated by a component, which 
has a specific shape roughly described by a broken power-law. The low-frequency part is 
mostly flat, while the power-law index a above the “break” frequency u br is variable between 
1 and 2. It is well established that the fractional root-mean-square (rms) variability in a 
source lightcurve decrease as a source evolves from LH state to HS state. Simultaneously, 
both u br and a increase. Although empirical shot-noise models were able to describe in 
general the observed PDS shape (Focke et al. 2005), the physical picture explaining the 
observed evolution during spectral transitions was missing. Moreover, shot-noise models 
were challenged by the linear absolute rms-flux relation (Uttley, McHardy & Vaughan 2005; 
Uttley 2004). This rms flux relation assumes that amplitudes and time-scales of shots are 
not independent, but are related in some way. 

Lyubarskii (1997) was the first to suggest a model for this time variability production 
in the accretion powered X-ray sources. He considered small amplitude local fluctuations 
in the accretion rate at each radius, caused by small amplitude variations in the viscosity, 
and then studied the effect of these fluctuations on the accretion rate at the inner disc edge. 
His linear calculations show that if the characteristic time-scale of the viscosity variations 
is everywhere comparable to the viscous (inflow) time-scale, and if the amplitude of the 
variations is independent of radius, then the power spectrum of luminosity fluctuations is 
a power-law \/u. If the amplitude of the variations increases with radius, the slope of the 
power spectrum of the luminosity variations is steeper than 1. Lyubarskii pointed out that 
he had no physical model for the cause of such fluctuations. Uttley, McHardy & Vaughan 
(2005) pointed out that rms-flux relation is naturally explained in the framework developed 
by Lyubarskii. 

TSA07 formulated and solved the problem of local driving perturbation diffusion in a 
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disk-like configuration. The problem of the diffusive propagation of the space distributed 
high-frequency perturbations is formulated as a problem in terms of the diffusion equation for 
the surface density perturbations. This equation is combined with the appropriate boundary 
conditions. The formulation of this problem and its solution are general and classical. The 
parameters of the resulting PDS, diffusion time scale of the diffusion propagation of the 
local perturbations t 0 and the power-law index of the viscosity distribution over the disk-like 
configuration a, are essential parameters of diffusion in a given bounded configuration. In 
TSA07 we call our model for the PDS of Green’s function of the bounded configuration as 
a white-red noise (WRN) and we adopt this name throughout this paper. 

The problem formulation was similar to the Lyubarskii’s scheme. However, the method 
of solution is different. Lyubarskii treated the factorization of the driving term (i.e. separat- 
ing it into two parts each depending on time and radius only) by linearizing the system. In 
TSA07 the analytical solution is obtained for the case of the factorized driving sources. Then, 
using the mean value theorem we showed that the general solution is simply a convolution of 
the response (Green’s function) signal of the configuration and the mean driving signal in the 
configuration. Thus the resulting power spectrum of the X-ray signal, as a convolution, is a 
product of the power spectrum related to the disk- like configuration response (the Green’s 
function) and that related to the perturbation sources (sources of driving oscillations). The 
Green’s function PDS is a white-red noise power spectrum (WRN). Specifically, the low fre- 
quency (LF) asymptotic form of the WRN PDS, when the frequency is less than the inverse 
of diffusion timescale in the disk t^ 1 , is characterized by a flat shoulder (white noise). In 
other words, the LF white noise shoulder is insensitive to the source distribution and to 
the viscosity law in the disk as a function of radius. The high frequency (HF) asymptotic 
form of WRN is a power law u~ a with index a, which is determined by the viscosity and 
perturbation source distribution over the disk. When the viscosity linearly increases with 
radius and the perturbation sources distribution is quasi- uniform, the index is a = 3/2. The 
basis of the presented power spectrum formation scenario is that the timing signal of the 
WRN PDS shape is a result of diffusive propagation of driving perturbations in the bounded 
configuration' (disk or Compton cloud) in the same way as X-ray photon spectrum is a re- 
sult of the photon diffusion (namely, upscattering of seed photons) in the same bounded 
configuration. 

The driving oscillation amplitude is assumed to be a smooth function of the radius. 
TSA07 suggested that driving fluctuations in the configuration can be introduced by g-mode 
driving oscillations at any given disk annulus. The local g-mode driving fluctuations, pro- 
duced possibly by local Rayleigh- Taylor local instabilities are high-frequency damped quasi- 
periodic oscillations (QPOs) which frequencies are related to the local Keplerian frequencies. 
As we mentioned above TSA07 formulated and solved a problem of the diffusive propagation 
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of the space distributed high-frequency perturbations in the bounded configuration. 

Our diffusion model for PDS is a product of WRN PDS and the driving source PDS 
(Lorentzian). 

The WRN PDS is a power spectrum of the solution of the initial value (Cauchy) problem 
which is a linear superposition of exponential shots [see Wood et al. (2001)]. For example, if 
the driving perturbations are distributed according to the first eigen-function of the diffusion 
operator then the bounded medium works as a filter producing just one exponential shot as 
a result of the diffusive propagation of eigen-function distribution of the seed perturbations. 
In the general case the resulting signal is a linear superposition of exponential shots which 
are related to the appropriate eigen-functions. Furthermore, TSA07 demonstrate that the 
observed rms-flux relations [e.g. Uttley, McHardy & Vaughan (2005)] is naturally explained 
by our diffusion model. In the framework of the linear diffusion theory the emergent pertur- 
bations are always linearly related to the driving source perturbations through a convolution 
of the disk-response (Green’s) function and source distribution. 

An important question is what our diffusion model predicts for relative contributions of 
the WRN PDS and the driving oscillation PDS in the resulting PDS and for a dependence 
of the integrated PDS power on the driving oscillation frequency. The next question is how 
this model dependence of the integrated power vs the driving oscillation frequency is related 
to the observed dependence of that. The answers to these questions are the points of the 
presented study. 

In §2 we refer to details of Cyg X-l observations with RXTE. In §3 we outline the main 
features of the diffusion model and related formulas. In §4 we show how the model integrated 
power vs the driving oscillation frequency fits X-ray data from Cyg X-l. In §5 we present 
the inferred correlation of Reynolds number with the driving oscillation frequency and the 
spectral state (photon index). Application of the paper results to the observed index-QPO 
frequency correlations is discussed in §6. Conclusions follow in §7. 


2. Observations 

For our analysis we used Cyg X-l data from the Proportional Counter Array (PC A) and 
All-Sky Monitor (ASM) onboard RXTE [Swank (1999)]. The data are available through the 
GSFC public archive 1 . In this Paper we present the analysis of a representative subset of 
RXTE observations of Cyg X-l. A reader can find the details of data reduction and analysis 


fiittp:/ /heasarc. gsfc.nasa.gov 


-5- 


in Shaposhnikov & Titarchuk (2006) and TSA07. We chose approximately 200 observations 
to cover the complete dynamical range of the source evolution from low/hard to high/soft 
state. For the presented analysis we refit PDSs with our new model and we used the results 
of our previous spectral analysis for photon index T. 

To fit a PDS we used a sum of our perturbation diffusion model and one or two 
Lorentzians to account for the Low Frequency Quasi-Periodic Oscillations (LFQPOs). For 
higher photon inicies, when the source is close to high/soft state, the contribution of the 
accretion disk variability component sometimes becomes significant. It is observed as an 
additional power law at the lower frequencies (see TSA07 for details). We fit this component 
with simple power law, when it is needed. 


3. The main features of the model 


The resulting variability of X-ray signal is determined by the fluctuations of the lumi- 
nosity A L x (t). We assume that the mass accretion rate variations AM(0, t) is converted 
with efficiency e e ff into the variations of the X-ray luminosity, i.e. A L x (t) = £ e //AM(0,f). 
TSA07 show that the fluctuations of the resulting X-ray oscillation signal A L x (t) due to the 
diffusion of the driving perturbations is 


A L x 



<p(t')Y (t — 


( 1 ) 


i.e. a convolution of the response (the Green’s function) of the disk-like configuration Y(t) 
(WRN) and the source variability function ip(t). The resulting power spectrum is 

||F x M|| 2 = ||F,(n;)|| 2 ||F y (a;)|| 2 (2) 


where F x (u), F^iu), Fy(uj) are Fourier transforms of A L x (t), <p(t), Y(t) respectively. Using 
the total power of the driving oscillations Pd r one can present the driving oscillation PDS as 
(see Eqs. 22 and B5 in TSA07) 


\\F»\\l 


f 1 drP ( 'dr / ipl 0 

{V - ^r) 2 + (f *./2)2 


( 3 ) 


where Vd r is the driving oscillation frequency, T*- is a full width of half maximum (FWHM) 
of the Lorentzian and a constant a varies in the range between 1 and 2 depending on the 


ratio of 2u dr ./F dr : 


„ 2 arct an (2u dr /F dr ) 

a = 1 H 


( 4 ) 
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For example a = 1 and a = 2 when 2v dr /Y dr <C 1 and 2v dr /Y dr >• 1 respectively. 


In Fig. 1 we show a typical example of the model fit to the data using Eq. (2). The 
parameters of PDS continuum for white-red noise component (WRN) ||Fy(u;)|| 2 (see for 
details TSA07) are the diffusion time scale to, index of the power-law viscosity distribution 
0 and for the driving oscillation component ||F^(u>)|| 2 they are u dr and Y dr (see Eq. 3). 

The power-law index of the viscosity distribution 0 is related to the power-law index of 
the red noise in the WRN PDS (see TSA07): 


and 


3 . 16-50 

Q = 2 - 6 = WT) 


for 0 > 0 


a = 2 for 0 < 0. 


( 5 ) 


The model predicted integrated power is (see TSA07) 

Px ’ dlff ~ 270/^0(3 +1/4Q)DC (6) 

where Q = Y dr /u dr is a quality factor for the driving signal and D is a factor of order of unity. 
Equation (6) was derived use the mean value theorem for the integral of the product of two 
functions. In TSA07 we assumed that a constant C related to the mean value of ||Fp(i/)|| 2 
over the frequency integration range is about a few (see Appendix B2 in TSA07 for details). 
Here we specify and obtain C— constant when we compare the model dependence P x ,diff on 
v dr with the observable P x on u dr . In order to make this comparison one should determine 
the best-fit parameter to and P dr as functions of u dr . Note that to derive Eq. (6) we also use 
the fact that the WRN PDS ||Fy(u;)|| 2 is normalized to l/(Dto) where DZ 1. (see Eq. B16 
in TSA07). 


Now we follow the method suggested in TSA07 to infer P dr vs u dr from the observations. 
Namely, given the fact that the driving PDS is a constant at frequencies u <C v dr we have 




fjrPdrlia^) 

d + M 2 ' 


Because for any power spectrum ||F(a;)|| 2 


( 7 ) 


||F(u;)|| 2 dj = ||F(27rzz)|| 2 27rdz/ = \\F(u)\\ldu 


we obtain that (compare with Eq. 2) 


||F r ( i ,)|| 2 =,(27r)- 1 ||F,»|| 2 ||F y (zz)|| 2 . 


( 8 ) 
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Thus a combination of Eqs. (7) and (8) leads us to determination of the integrated 
power of the driving oscillations 


Pdr = 


+ (r*-/2) 2 ) 


ll^(0)||’ 


2t r 2 a^- 


+ (f dr /2) 2 \\F X (0)\? 


dr 


dr 


ll^(o))ll 5 


( 9 ) 


We remind a reader that the values of u dr and f*. are the best-fit PDS parameters, ||F X (0)| | 2 
is the observed PDS value at v = 0 and ||Fy(0))|| 2 is a value of the normalized WRN PDS 
at v = 0. As we mention above the integral of the normalized WRN PDS \\Fy(u))\\ 2 over v 
is 1 /(Di 0 ). 


4. The integrated power vs driving oscillation frequency and photon index 


In Figure 2 we present the observed correlation of low frequency QPO centroid v L with 
the driving QPO frequency v dr (upper panel) and photon index P with v dr (lower panel). 
These correlations imply that iq, along with v dr increase when the source becomes softer. In 
other words, the emission area [Compton cloud (CC)] contracts when the source evolves to 
the soft states. 


In Figure 3 we also see this effect of CC contraction as anticorrelation of the CC diffusion 

time scale t 0 with v dr . The inferred dependence of t 0 vs u dr can be fitted by the power law 

j. „ ..-2.13i0.14 
to oc u dr 

We infer the integrated power of the driving oscillations P dr vs u dr (see Eq. 9) and 
then we obtain the model integrated power P x ,diff vs u dr (see Eq. 6). In Figure 4 we show 
that the dependence of P dr vs v dr can be fitted by power law P dr oc ^ 1 ' 8±0 ' 16 . Namely 
the driving oscillation power P dr decreases when the source (Cyg X-l) goes to softer states. 
Presumably the decay of P dr with v dr is also related to the contraction of Compton cloud. 
The driving oscillations can result from the Rayleigh-Taylor (RT) local instability [see e.g. 
Chandrasekhar (1961) and Titarchuk (2003)] ]. The decay of P dr can be considered as a 
cumulative effect of the local Rayleigh-Taylor (RT) instability when the effective area of a 
given configuration (CC) undergoing RT oscillations contracts. 

In Figure 5 we present a comparison of the observable PDS integrated power P x (black 
filled circle) with the model predicted P x , d iff (crosses) (see Eq. 6). One can see that the 
dependence P x ,diff on u dr is similar to the observable correlation P x vs v dr . 


Note that we obtain the factor C ~ 4 (see Eq. 6) by shifting a set of the values of P x , dl / / 
along Y-axis to fall on top of P x values. The power-law P x oc ^r°' 48± °' 03 fits the dependence 
of the theoretical and observable integrated powers vs the driving oscillation frequency. 
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5. Reynolds number of the accretion flow in Compton cloud configuration 

Using the best-fit parameters of the PDS we can also infer the evolution of the physical 
parameters of the source such as Reynolds number of the accretion flow Re, with the change 
of photon index T. In fact, TSA07 relate t 0 with Re and a magneto-acoustic QPO frequency 


z'M a as 


to = -■ 


VmaR o 


R o 


Re 


3 (4 - ip ) 2 [ v(Ro) J \VmaJ 3 (4 - ip ) 2 a MA v MA 
where clma is a numerical coefficient. Formula (10) leads to the equation 


( 10 ) 


Re = 27r (^f)| - 4 — (^o) ( U > 

that allows us to infer a value of Re using the best-fit model parameters to and the QPO low 
frequency Vl which presumably equals to vma- Ultimately we can find the evolution of Re 
with the photon index P given that vl, to and the viscosity index ip evolve with T (see Figs. 
2-3 and 6). 


In Figure 7 we present the inferred Reynolds number as a function of the photon index 
T. We use Eq. (11) where we set a AA [ = 2tx (see details of this assumption in TSA07) 
and the observable correlations of v L and t 0 with F (see Figs. 2 and 3). One can see that 
Re— number steadily increases from 10 to 70 when the source evolves from low/hard state 
to high/soft state. In contrast, TSA07 found that Re ~ 8 ± 2.5. We note, however, that in 
TSA07 only WRN model was used without accounting for driving oscillation distribution, 
which significantly affects the resulting value for ip. They also used a limited set of data. 


Note the observed behavior of the Re-number vs F and mass accretion rate was predicted 
by Titarchuk, Lapidus & Muslimov (1998), hereafter TLM98, where they formulated a 
transition layer model (TLM) and studied its consequences for observations. It is important 
to emphasize that the Re— number along with the photon index F, the low frequency QPO 
v L and the driving frequency v dr can be considered as characteristics of the spectral state. 
All of them correlate with each other. 


6. Photon index-QPO frequency correlation 

TLM98 showed that the outer (adjustment) radius of the transition layer R out mea- 
sured in the dimensionless units with respect to Schwarzschild radius Rs = 2GM/c 2 , r out = 
Rout/Rs, anticorrelates with Re— number or photon index F (spectral state) only. Thus vl 
(or vj[j A ) as a ratio ~ VjU A /R 0 ut should correlate with T (or Re— number) where values vl 
related to the same F for different sources should be inversely proportional to a mass of 
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the central object (black hole or neutron star). Note that a plasma velocity Vma is also a 
function of T only. The comparison of the observed index-QPO frequency correlations for 
two different sources (with two different masses) should lead to determination of their rela- 
tive masses with respect each other. This is the main idea behind the method of weighing 
black holes (TLM98) recently applied for BH mass determination in a number of Galactic 
and extragalactic sources [see Titarchuk, & Fiorito (2004); Fiorito & Titarchuk (2004); 
Dewangan, Titarchuk &; Griffiths (2006); Strohmayer et al. (2007) and Shaposhnikov & 
Titarchuk (2007)]. 


7. Conclusions 

We explain the decay of the emergent time variability of X-ray emission in compact 
sources when these sources evolve from low/hard to high/soft states. We find that the 
resulting power P x from Cyg X-l decays with the driving oscillation frequency v dr as P x oc 

We show that the reciprocal of the diffusion time scale of the perturbation 1 , the low 
frequency QPO u L , the driving oscillation frequency u dr , inferred by fitting the Cyg X-l 
PDSs with our diffusion model, increase when the source evolves from low/hard state to 
high/soft state. This behavior of the PDS characteristics implies that the Compton cloud 
contracts towards softer spectral states. The driving oscillations are probably caused by the 
local Rayleigh- Taylor instability which cumulative P dr decreases when the effective area of 
the configuration producing the RT oscillations contracts. The decay in driving power leads 
to the decay in the total observed variability power from the source. Using the fact that 1 , 
Ul, u dr increase with T and P dr , P x decrease with F we conclude, as a result of our analysis, 
that the Compton Corona shrinks when Cyg X-l goes from low/hard state to high/soft state. 

Our extensive data analysis of the power spectra from Cyg X-l with an application 
of the method of Re— number determination developed in TSA07 indicates that Re related 
to Compton cloud configuration increases from values about 10 in low/hard state to that 
about 70 in high/soft state. We confirm the predictions by TLM98 that Re— number should 
increase with index and QPO frequencies. Thus one can conclude that the observable index- 
QPO correlation is probably driven by the increase of Re— number when the source evolves 
from low/hard state to high/soft state. It is worth noting that inverse proportionality of the 
low-frequency QPO with respect BH mass in the index-QPO correlation leads to the method 
weighing BHs employing this index-QPO correlation. 
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Fig. 1. — A particular example of observable PDS. The PDS continuum is fitted by our 
diffusion PDS model which is a product of WRN PDS and the driving oscillation Lorentzian. 
We also use a simple Lorentzian to fit QPO features. 
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best-fit model parameter, diffusion time scale t 0 vs u dr . The dashed line is the 
law t 0 oc v~[ r 213±0U . 
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Fig. 4. — Inferred P<* r vs ^ r . The dashed line i 



best fit power law P^ r oc ^ dr L8±0 ' 16 . 




v dr , Hz 

Fig. 5. — Comparison of the model P Xtdl ff vs u dr (crosses) (using Eq. 6) with the observable 
P x vs v dr (black filled circle) The dashed line is the best fit power law P x oc z/7°' 48±0 ' 03 . 
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